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Abstract. We describe a double-peak structure of the coincidence spectrum of 
biphoton states in the process of spontaneous parametric down-conversion with a pump 
having the form of short pulses. The effect is shown to arise owing to the obligatory 
symmetry of bihoton wave functions, as wave functions describing states of two bozons 
obeying the Bose- Einstein statistics. Parameters of the peaks are found and conditions 
necessary for experimental observation of the effect are determined. 
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1. Introduction 



As known [1, 2], biphoton states generated in Spontaneous Parametric Down- Conversion 
(SPDC) can be highly entangled in continuous variables of signal and idler photons, such 
as frequencies or angular variables. An efficient way of finding the degree of entanglement 
from experimental data consists in measuring single-particle and coincidence angular 
distributions or spectra, finding their widths (correspondingly, Au;^*^ and Aa;^"^^ for 
spectral distributions) and using their ratio R = Au^^^ Aco^'^^ as the entanglement 
quantifier. As shown in Refs. [3, 4], for double- Gaussian distributions R = K, where K 
is the Schmidt parameter [5, 6]. Both parameters K and R can be used as entanglement 
quantifiers only for pure biphoton states. The approximate equality R ^ K assumes 
that both single-photon and coincidence distributions have smooth one-peak forms and 
can be approximated by Gaussian curves. In this paper we show that there are cases 
when both assumptions about purity of biphoton states and about Gaussian-like shape 
of the coincidence distributions are invalid. The key point of this consideration is related 
to the obligatory symmetry of biphoton wave functions with respect to transpositions 
of the photon variables. Importance of this feature of any two-particle states with 
identical (indistinguishable) particles was discussed in details in Rcf. [7] in application 
to systems with discrete variables (biphoton qutrits and ququarts). Here we consider 
the case of continuous variables (frequencies of photons) in the SPDC processes with 
coUinear phase-matching and with the pump having the form of a short pulse. 



1.1. Type-I coUinear phase-matching 

Let us consider first the case of the typc-I coUinear phase-matching. In this process 
the pump propagates in a nonlinear crystal as an extraordinary wave, and some of its 
photons decay for two photons propagating in the same direction (z-axis) as ordinary 
waves, e 0-1-0. Both emitted photons have the same polarization, e.g., the horizontal 
one, whereas the pump is vertically polarized. We assume also that the pump is a 
Gaussian pulse of a duration r. Under these conditions the polarization-spectral wave 
function of two emitted photons has the form [2] 

where ui and uj2 arc frequencies of emitted photons, ujq is the central frequency of the 
pump, is the angle between the optical axis of a crystal and the z-axis (direction of 
propagation of all photons), L is the length of a crystal in the propagation, z, direction. 



Two columns in Eq. (1), (^gj^ and (^0/2' polarization parts of the biphoton wave 
function with indices 1 and 2 referring two emitted photons. A{uJi,uj2,(j)) is the phase 
mismatch 

A{uJi,UJ2,(p) = kp{ui +UJ2,4>) - ki{ui) - k2{uJ2) 

= ne(a;i -Fa;2,0) ^ Uoiuji) — -no{uj2) — , (2) 
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where and no(a;) is the isotropic index of an ordinary wave (emitted photons), and 



is the refractive index of the pump 



ne(w,0) 



no(w)ne(w) 



[{no{uj) sin 0)2 + (ng 



, uo cos 



0)2]l/2 



(3) 



ne{(jj) is the refractive index of an extraordinary wave that would propagate along the 
optical axis. The functions no{(jj) and ne{(jj) are determined by the well known Sellmeier 
formulas [8]. 

Without any further detailing, it's clear from Eqs. (1) and (2) that in the case of the 
type-I phase-matching the polarization-spectral biphoton wave function is symmetric 
with respect to the transposition of particle variables 1^2 and does not need any 
additional symmetrization. Actually, symmetry of the wave function (1) is related to 
the fact that both emitted photons belong to the same type of waves, ordinary waves. 

Another general comment concerns the polarization degree of freedom of emitted 
photons. The polarization part of the wave function (1) is factorized for the product 
of two terms depending separately on polarization variables of photons 1 and 2. Also, 
it is factorized with respect to the part depending on frequency variables. For this 
reason, in the SPDC process with the type-I phase-matching, there is no polarization 
entanglement. Moreover, a pure biphoton state characterized by the wave function of 
Eq. (1) remains pure even if it is reduced with respect to the polarization variables. 
In other words, the bihoton state has only purely frequency entanglement, and for 
its investigation polarization part of the wave function simply can be ignored. This 
contrasts with the case of the type-II phase-matching considered in the next section. 




Figure 1. Angular dependence of the phase mismatch at central frequencies, 
A(cjo/2,cjo/2,(/)) . 



To specify a little bit further the most typical experimentally met conditions, let us 
assume that the SPDC process under consideration is degenerate with respect to central 
frequencies of spectra of emitted photons, ujf^ = ^2°'' = i^o/2- Such situation occurs 
if at these frequencies the phase mismatch turns zero, A{ll!o/2,uo/2, (p) = 0, which can 
be provided by an appropriate choice of the angle (p. Figure 1 shows the dependence 
of the phase mismatch at central frequencies on 0. Here and below all calculations are 
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made for a crystal LilOa and the central wavelength of the pump Aq = 0.3975/im. As 
can be found from Figure 1, the degeneracy condition A (a;o/2, Wo/2, 0) = is satisfied 
at = 00 = 42.904°. 

Figure 2 shows a typical spectral dependence of the mismatch A at = 0o 
but frequencies different from those providing exact phase-matching. In this and all 
further figures all calculated spectra are plotted in dependence on wavelengths rather 
than frequencies, which is more appropriate for comparison with existing and possible 
experiments. As seen from Figure 2, the mismatch is a monotonous function of Ai at 



a given value of A2 (or of a frequency ui at a given uj2)- For this reason, the mismatch 
turns zero at only one value of Ai (or uji = 27rc/Ai). As the zero mismatch corresponds 
to the maximal value of the sinc-f unction in Eq. (1), the coincidence spectrum has in 
this case a simple single-peak structure (Figure 3), in contrast to the spectrum arising 
in the case of the type-II phase-matching discussed below. 
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Figure 2. Mismatch (2) as a function of Ai at A2 = 0.9/im. 
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Figure 3. Coincidence spectrum, |\I>(Ai,A2)P (in arbitrary units), at A2 = const = 
0.794/im; the type-I collinear degenerate phase- matching, pulse duration and central 
wavelength of the pump r = 50/s and Ao = 0.795/xm, LilOa crystal of the length 
L = 1cm. 
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2. Type-II collinear phase-matching 

In the case of the type-II phase-matching, the pump and one of two emitted photons 
are extraordinary whereas the second emitted photon is of the ordinary waves, i.e., the 
SPDC process has the form e — )■ o + e. As previously, let us assume that the process 
is collinear and degenerate with respect to the central frequencies of all photons. Let, 
as previously, polarizations of extraordinary and ordinary photons be, correspondingly, 
vertical and horizontal. The phase mismatch of such process at central frequencies is 
given by 

A(^^)(0) = fcp(a;o/2,0) - K{ujom - fce(c^o/2» 

= n.K, 0)^ - (^) ^ - (f , 0) |. (4) 

For LilOs and Aq = 0.39975/im the function A*^^^)(0) is plotted in Figure 4. The exact- 




Figure 4. Angular dependence of the phase mismatch at central frequencies 
A(^^)(0) (4). 



phase-matching angle at which /S'^^^^cf)) = found from Figure 4 is 



4<"' 



Both this value and the curve of Figure 4 rather pronouncedly differ from 



68.845°. 
.0 - 42.904° 

and the curve of Figure 1 occurring in the case of the type-I phase-matching. 

At frequencies of emitted photons Ui and 002 differing from ijJo/2, the phase 
mismatch A of the type-II SPDC process is asymmetric with respect to transposition 
of photo variables (see Eqs. (6) and (7) below). The same is true for wave functions 
constructed from a given A. But as photons are bozons, their wave functions must 
be symmetric [9]. Symmetrization of wave the function is provided by summation of 
its expression with a given A and the same expression with transposed frequency and 
polarization variables of photon : 

{oOi + UJ2 



"if oc exp 




81n2 



smc 




smc 



L ^ 
-A21 



(5) 



where 



A 



12 



(Wi, UJ2) = kp (ui + UJ2, (po^^) - koM - ke (W2, 0, 
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= He [UJI + UJ2, (pQ ) no{UJi)— - ne[uj2,(po — , (6) 

A2l(Wl, UJ2) = kp (ui + LO2, (f)o^^) - ko{uj2) - ke (uji, 0^1^^^) 

= He [Ui + UJ2, (po ') no{UJ2)— - ne{uJi,(p'Q ')—. (7) 

Mismatches Ai2{uJi,U2) and A2i{uJi,U2) are equal to each other only at ui = 002- 
Otherwise they are different. In particular, in dependence on one of the frequencies 
(e.g., Ui) at a given value of another one {002) mismatches A12 and A21 turn zero at 
different values of Ui. This is illustrated by Figure 5 where, as usual, A12 and A21 
are considered as functions of wavelengths rather than frequencies. 5). Appearance 
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Figure 5. Mismatches A12 (6) and A21 (7) as functions of Ai at A2 = 0.9fim for LilOa 
and Ao — 0.795/xm . 

of two different values of Ai (or, equivalently, ui) indicates a possibility of observing 
two different peaks in the coincidence spectra. But, actually, one can observe either 
single-peak or double-peak coincidence spectra depending on conditions of observation. 

Indeed, the usual way of observing coincidence spectra consists in splitting the 
SPDC biphoton beam by a nonselective beam splitter for two channels A and B with 
identical beams in each channel. Detectors (photon counters) are installed in each 
channel, and signals from each detector are sent to the computer to select only coinciding 
events, i.e., counts of photons of the same rather than different SPDC pairs. Parameters 
of one detector (e.g., in the channel B) are kept constant whereas parameters of the 
other detector (A) are scanned to reconstruct is such way the coincidence distribution 
of photons (number of counts) over the varying parameters. For finding spectral 
distributions detectors must be spectrally sensitive, i.e., having frequency filters installed 
in front of photon counters. If in addition detectors are provided with polarizers in 
front of photon counters, detectors are polarization-selective. In this case one can see 
separately only one of two coincidence peaks determined by the wave function (5). For 
example, if the detector in the channel B counts only vertically polarized photons, and 
we denote their frequency as U2, these photons belong to the type of extraordinary waves 
and are described by the first term in the sum of Eq. (5). Under these conditions the 
detector in the channel A will count only horizomtally polarized (ordinary) photons, and 
its scanning over frequency Ui (or wavelength Ai) will reproduce the single coincidence 
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spectral peak around the frequency (wavelength) where A12 = 0. Oppositely, if the 
detector in the channel B is provided with the horizontal polarizer, the second detector 
(in the channel A) will register only vertically polarized photons, to reproduce the only 
coincidence peak described by the second term in Eq. (5) and located near the frequency 
(wavelength) where A21 = 0. To see both coincidence peaks together, one has to remove 
polarizers and to make detectors frequency-selective but non-sensitive to polarizations 
of photons. 

The last case of measurements by only frequency-selective detectors corresponds in 
theory to construction of the density matrix from the wave function (5) and reduction 
(tracing) of this density matrix over the polarization variables. The arising frequency 
density matrix is given by 



pioji, U2; oj'i, 002) oc exp 
X I sine ^Ai2(a;i,a;2) 



81n2 



smc 



-Ai2{uj[,UJ2) 



-|-smc —A2i{oJi,uj2) smc 
2 



2^2i(^,^) 



(8) 



This density matrix describes a mixed rather than pure state with continuous variables. 
Even definition of the degree of entanglement for such states is a problem having no 
clear and recognized solution. But measurable distributions of probabilities are rather 
easily found because they are directly determined by diagonal elements of the density 
matrix (8). The two- frequency probability density is given by 



dw 



diij\duj2 



= p{ui,(jj2.,0Ju0J2) OC exp \ -{Ui +U2- OJoY 



41n2 



X < smc 



'-Ai2{i^l,i^2) 



+ sine 



i2l(<^l,C^2j 



(9) 



Conditional probability density for one photon to have a frequencies uii under the 



I u}=const. 1 



is 



condition that the other photon has some given frequency a;2, dw^'^^dui 
determined by the same equation (9). Alternatively, this equation can be rewritten in 
terms of wavelengths: 



dXi 



X < smc 



exp 



L ^ 
- A12 
2 



1 

2ttc 



1 

2nc 



1 



41n2 



Ai A2 



-|- sine 



L ^ 
- A21 
2 



2-RC 2-KC 



Ai A2 



(10) 



Of course, the theoretical concept of the conditional probability density is identical 
to the experimentally measurable coincidence distribution. Two typical examples of the 
double-peak coincidence spectra are shown in Figures 6a and 7. These two pictures 
differ from each other by different pulse durations of the pump. As seen, in the case of 
longer pulses (Figure 7) the double-peak structure is spectrally compressed compared 
to the case of longer pulses (Figure 6a). The picture of Figure 66 shows how two peaks 
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of the double-peak structure merge into a single peak when wavelengths (frequencies) 
of two emitted photons become equal to each other. 




0.79 0.795 0.8 0.805 0.81 0.815 0.793 0.795 0.797 



Figure 6. Coincidence spectra, dw^'^'>/dXi (in arbitrary units) at A2 = const. = 
(a) 0.79/zm and (&) 0.795/xin; the type-II collinear degenerate phase-matching, pulse 
duration and central wavelength of the pummp r = 50/s and Aq = 0.795/im, LilOa 
crystal of the length L = 1cm. 
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Figure 7. The same as in Figure 7(a), but at t = 186 fs and A2 = 0.7925 /im. 



3. Heights and widths of peaks and temporal walk-off 

Practically always both coincidence and single-particle SPDC spectra are narrow 
compared to the pump frequency uo- For this reason wave vectors in the expressions (6) 
and (7) for phase mismatches can be expanded in powers of deviations of frequencies ui 
and U2 from their central meanings Ui = Ui — Uq/2 and 1/2 = u)2 — ujq/2: 

kp{ui + U2) ^ kp{uo) + (z/i + U2)kp = kp{uo) + '^^-^ 

V\ 2 

ke{(^l,2) ~ fce(Wo/2) + Ul,2K = ^e(c^o/2) + -7^, 

^0(^1,2) ^ ko{uJo/2) + Vi^2k'o = ke{uJo/2) + -7^, (11) 
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where y^P'°'^') are the group velocities of the pump and emitted ordinary and 
extraordinary photons 



rip) 



dkp{oj) 



du 



-1 



(o,e) I dko,e{uj) 



W=£J0' 



duj 



(12) 



a;=aJo/2 y 



In terms of 1/1^2 the wave function of Eq. (5) takes the form 



"if oc exp 



81n2 



smc 



2c 



+ 



smc 



— {AoU2 + AeUi) 

Jc 



where A^ and Ag are the temporal-walk-off constants 

f I 1 \ . f I 1 



(13) 



(14) 



For LilOa and Aq = 0.3975/xm the group velocities v^'"'^^ and the walk-off constants 
Ao^g are easily calculated to be given by: wj^') = 0.4986c, v^f^ = 0.522c, v^'^) = 0.5628c, 
^0 = 0.09, and Ae = 0.2287. 

Let us assume now that that the pulse duration of the pump is so short that 
the Gaussian function in Eq. (13) is much wider than both sinc-functions in their 
dependence on i/i, i.e., that 

T <C fti 3ps. (15) 

c 

Then, positions and widths of the coincidence spectral peaks arc determined mainly 
by the narrow sinc-functions. In dependence on Ui at a given U2, peaks related to 
mismatches A12 and A21 are located at and z^j^^' and have widths Az/I^^"* and 
Avi^^^ given by 



(12) 



(12) 



f= n (21) 

-^^2-7- = -2.54z/2, i^i 

Aq 



-^2-7- = -0.39i/2, 



LA 



LA 



(16) 
(17) 



In terms of wavelengths the same relations take the form 



(12) 



2Ao = -(A2 - 2Ao)-^ = -2.54(A2 - 2Ao), 



A^ ^ — 2Ao — — (A2 
AA^) 



2Ar 



\2 

AX (21) 



3^2 



-0.39(A2-2Ao), 



(18) 
(19) 



LAg LAg 

These equations show that the peaks of the coincidence spectrum dw/dcui are located 
always at the opposite side of the central frequency iOo/2 or wavelength 2Ao compared 
to the frequency U2 (wavelength A2) of photons registered by the detector with fixed 
parameters (in the channel B). The peak located closer to the central frequency 
(wavelength) is described by the term with the mismatch A21 in Eq. (5) and is 
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significantly narrower than the second wider peak, located almost ten times further from 
the central frequency (wavelength) and related to the term with A12 in Eq. (5). Both 
the double-peak structure and its described features arise owing to the well pronounced 
difference between the walk-off constants and Ae, and this difference originates 
from the fact that two emitted photons belong to different types of waves, ordinary 
ad extraordinary, with different group velocities. If the constants A^ and A^ would be 
equal, two peaks in the coincidence spectrum would merge into one, and this is the 
situation occurring in the type-I SPDC process, in which both emitted photons are 
ordinary waves. 

As for heights of the peaks in the coincidence spectrum, as long as we do not 
calculate normalization factors in wave functions and density matrices, absolute values 
of the peak heights cannot be found. But in arbitrary units, in the approximation of 
short pump pulses (15), the peak heights are determined by the Gausssian function in 
Eqs. (5), (8), (9), (13) in which ui is substituted either hy Uo/2 + ' or by Ci;o/2 H- z/j : 



^ / max V 



12) \2 I [A,.-Anrv^T- 



>. ^ / max V 



exp 




The ratio of peak heights equals to 

idw/duj^)(^^ i {Al - A^:){A, - A^f vIt^ \ 

If this ratio is large, the farthest peak becomes practically invisible. On the other hand, 
the ratio (22) can be minimized by making the fixed frequency 1^2 sufficiently small. 
The question that arises in such case is how small is a reasonable value of 11^21? The 
answer is evident: the minimal value of |z/2| is that at which the distance between peaks 
becomes on the order of width of the farthest peak, i.e., |z/2|min ~ X ^2*1^2 ■ Substitution 
of this value into the right-hand side of Eq. (22) gives the minimal achievable value of 
the peak-height ratio: 



_{dw / dLJi)mll 



max J jjjjjj 



= exp —————— -— . (23) 



Al{A, + A,) 4L2 1n2 

Under the condition (15) the expression under the symbol of exponent in this equation 
is small and, hence, the peak-height ratio is close to unit. This means that, indeed, 
in the femtosecond region of pulse durations of the pump the double-peak structure of 
the coincidence spectrum is rather well pronounced in a sufficiently large region of fixed 
frequencies |i/2| exceeding |f2|min and is accessible for experimental observation. 



4. Conclusion 



Thus, the effect described above consists in appearance of a double-peak structure in the 
coincidence spectrum of the biphoton state arising in the SPDC process with a coUinear 
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type-II phase-matching, short duration of pump pulses and degeneracy with respect 
to central frequencies of emitted photons. The effect is shown to be related directly 
to the symmetry of the biphoton wave function with respect to the transposition of 
photon variables. Though often not noticed or ignored, this feature of two-bozon wave 
functions is shown to be responsible for a clearly observable and not negligible effect of 
peak doubling. The conditions under which the effect can be observed are determined. 
The main of them consists in the requirement that the pump pulses have to be short 
enough, of a femtosecond duration. Positions, heights and widths of the peaks are found 
and expressed in terms of the temporal-walk-off parameters of the SPDC process. 

5. Acknowledgement 

The work is supported partially by the RFBR grants 08-02-01404, 10-02-90036, and 
11-02-01043 

References 

[1] M.V. Fcdorov, M.A. Efermov, P.A. Volkov, E.V. Moreva, S.S. Straupe and S.P. Kulik, Phys. Rev. 
A, 77, 032336 (2008) 

[2] Yu.M. Mikhailova, P.A. Volkov and M.V. Fedorov, Phys. Rev. A, 78, 062327 (2008) 
[3] M.V. Fedorov, M.A. Efremov, A.E. Kazakov, K.W. Chan, C.K. Law, and J.H. Eberly, Phys. Rev. 
A, 69, 052117 (2004) 

[4] Fedorov M V, Efermov M A, Volkov P A and Eberly J H, 2006, J. Phys. B, 9 S467 
[5] R. Grobe, K. Rzazewski, and J. H. Eberly, J. Phys. B, 27, L503 (1994) 
[6] A. Ekert and P.L. Knight, Am. J. Phys., 63, 415 (1995) 

[7] M.V. Fedorov, P.A. Volkov, Yu.M. Mikhailova, S.S. Straupe, S.P. Kulik, arXiv:1009.2744v2 [quant- 
ph] 27 Sep 2010 

[8] W. Sellmcicr, Ann. Phys. Chem., 143, 271 (1871) 

[9] A.I. Akhiezer and V.B. Berestetskii, Quantum Electrodynamics, Interscience Publishers 1965, a 
division of John Wiley & Sons, New York-London-Sydney [Chapter I, § 7, Eq. (7.3)] 



